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We discuss the relevance of a dedicated measurement of exclusive production of a pair 
of neutral pions in a hard 7*7 scattering at small momentum transfer. In this case, 
the virtuality of one photon provides us with a hard scale in the process, enabling us to 
perform a QCD calculation of this reaction rate using the concept of Transition Distribution 
Amplitudes (TDA) . Those are related by sum rules to the pion axial form factor F\ and, as 
a direct consequence, a cross-section measurement of this process at intense beam electron- 
positron colliders such as CLEO, KEK-B and PEP-II, or Super-B would provide us with 
a unique measurement of the neutral pion axial form factor F\ at small scale. 

1 Introduction 

In a series of papers pQ [31 El SI El El El [5] , we have advocated that factorisation theorems [H] 
for exclusive processes may be extended to the case of other reactions such as (Mj stands for a 
meson and Bi for a baryon) B\ B2 — > 7* 7, B\ B2 — > 7* Mi j T B\ — > B27, It —> B?.M\ or 
7^ 7 — > M1M2, in the kinematical regime where the off-shell photon is highly virtual (Q 2 of the 
order of the energy squared of the reaction) but the momentum transfer t is small. This enlarges 
the successful description of deep-exclusive 77 reactions in terms of distribution amplitudes [10] 
and/or generalised distribution amplitudes |llj on the one side and perturbatively calculable 
coefficient functions describing hard scattering at the partonic level on the other side. 

Intense beam electron colliders, such as B factories, are ideal places to study such reactions 

as 

7i7^P ±7r=F , ll 7^ 7r±7I " =F , 7l7 ~+ t^ 71 " , 

in the near forward region and for large virtual photon invariant mass Q. Recently BaBar 
reported a new measurement of the reaction 7*7 — > ir° up to photon virtualities squared of 
40 GeV 2 [T2J. In the latter study, the reaction 7*7 — * tt°tt q was investigated in the /2(1270) 
and /o(980) resonance region as a potential background for the study of the ir° transition form 
factor. This low-W 2 ^ kinematical region should be analysed in the framework of generalised two- 
meson distribution amplitudes |llj and in particular should solve the much discussed problem 
of its phase structure around the /o mass [13] which is of crucial importance for the ability to 
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detect Pomeron-Odderon interference effects in high energy electro-production of meson pairs 

We want here to emphasise another kinematical region, namely the small-i large- region 
at moderate Q 2 (2 GeV 2 and more) which may provide us with unique information on the tt° 
axial form factor at small scale which so far has never been experimentally measured. It has 
been argued that a new duality [TS] relates these two factorisation regimes. 

In principle, another possibility to study this quantity would be the crossed channel, that 
is DVCS on a virtual neutral pion along the lines exposed in Ref. [TB] for ir + . 



2 Pion-pair production in the TDA regime 

Let us recall the main ingredients of the analyses developed in [TJ H] focusing on the neutral pion 
case. With the kinematics described in Fig. [I] we define the axial 7 — > tt transition distribution 
amplitude (TDA) A(x, £,t) as the Fourier transform of matrix element (7r°(p n .)| |7(p 7 )) 
where Oa — tp(^f)[^, §] 7' 1 7 5 V'(f )■ The Wilson line [=f, §] ensures the QCD-gauge invariance 
for non-local operators and equals unity in a light-like (axial) gauge. We do not write the 
electromagnetic Wilson line, since we choose an electromagnetic axial gauge for the photon. 
We then factorise the amplitude of the process 7^7 — > 7r 7r° as 

£ / dxdz$l{z)Ml{z,x,t) q V , (1) 
, J fa 

q—u,a 

with a hard amplitude M?(z,x,£), ^(z) the distribution amplitude (DA) for q quark content 
of the tt meson with momentum p' n and (x, £, t) the axial 7 — > n TDA for the quark q. 




Figure 1: The factorised amplitude for 7*7 — > 7r°7r° at small transfer momentum. 

The variable z is as usual the light-cone momentum fraction carried by the quark entering 
the pion with momentum p' v , x + £ (resp. x — £) is the corresponding one for the quark leaving 



(resp. entering) the TDA. The skewness variable £ describes the loss of light-cone momentum 
of the incident photon and is connected to the Bjorken variable Xb- 

Contrarily to the case of generalised parton distributions (GPD) where the forward limit is 
related to the conventional parton distributions measured in the deep inelastic scattering (DIS), 
there is no such interesting constraints for the TDAs. The constraints we have here are sum 
rules obtained by taking the local limit of the corresponding operators and soft limits when the 
momentum of the meson in the TDA vanishes. 

Let us consider in more detail the 7 — > ir° axial TDAs which is defined by (P = p *~£ p "' ; 
A = Pit -p 7 ): 
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< ~ ) — ;- ri 5 qUh) = j^^(e-&)P tl A* (x,t,t) (2) 
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A sum rule may be derived for this photon to meson TDA by integrating on x both side of 
Eq. ([2| and we get 



e 



-(e ■ A)P» f dx Af(x,£,t) = (7r°|q(0) 7 ' 1 7 5 g (0)| 7 ), 

TT J — I 



(3) 



The latter matrix element of a local quark-anti-quark operator is directly related to the quark 
q contribution F\ to the axial form factor of the ir° meson. Similarly, we have in the vector 
charged pion case 0]: 

dxV^(x,U)^—F^(t), (4) 
1 m n 

with Fy" = 0.017 ± 0.008 Q2]. 

This sum rule constrains possible parametrisations of the TDAs. Note, in particular, the 
^-independence of the right hand side of the relation. 

3 Models and cross section evaluation 
3.1 Amplitude 

Let us thus consider the 7r°7r° production case when the ir° with momentum p' n flics in the 
direction of the virtual photon and the other ir° emerges from the TDA. For dcfinitcncss, we 

choose, in the CMS of the meson pair, p = 2 fi+£)w' (1 ? 0; 0,-1) and n — 2(02+^2") (■*-> 1) 
and we express the momenta trough a Sudakov decomposition (with = j^t and neglecting 
the pion mass): 

qi , W 2 r>2 a 2 

Pl = (l + 0p, 1= i + ™ n- q2 * w2 (I + OP, P w = (l-0P- I f^n + Ar. (5) 

We can see that £ is determined by the external kinematics through £ ~ ep^W 2 — similarly 
o 2 

to xb = q2+ W 2 to which it is linked via the simple relation £ ~ 2-x B ■ 
The hard amplitude amplitude in Eq. ([TJ thus reads : 

„, 9 , e . 4tt 2 a em a s C F Q q 1/1 1 \ 

M«(z,x,£) = — — q -—[ — - + —7 e-A, 6 

NcQ zz\x — t;+ie x + ^ — ieJ 



where Q u — 2/3, Qd = — 1/3 and with z = 1 — z. Note that the factor f„ in the ir DA &%{z) 
cancels with the one from the TDA definition and docs not appear in Eq. Q. Now, if we 
choose the asymptotic form for the neutral pion DA, <&^ (z) = — <5>^ (z) = *^z(l — z), the 
z-integration is readily carried out and after separating the real and imaginary parts of the 
amplitude, the ^-integration gives: 



A?(x,(,t)= (7) 



+ ie x + £ — ie 

7| 2^ IQ? I + A i (6&*)0°g(^|J -«0+ 

— — + A£ (-£, f , t)(log + itt) 



q—u,d 



V2 

q—u 



The scaling law for the amplitude is 

M™ A (Q 2 ,Z,t)~°^El, (8) 
up to logarithmic corrections due to the anomalous dimension of the TDA and the running of 



3.2 Remarks on available models 

Lacking any non-perturbative calculations of matrix element defining TDAs, we have initially 
built a toy model [I] based on double distributions [TB] to get estimates for the cross sections, 
to be compared with experimental data. In 0] , we compared the rate obtained with this model 
with the one from the model built in [19]. Subsequently, a model based on quark spectral 
representation was developed in [20] , another based on N JL model was studied in [2H H2] and 
lastly the ir — ► 7 TDAs were studied in a non-local chiral quark model [33]. All the models 
(see e.g. [23]) used so far for the pion GPDs could be extended to the construction of ir — ► 7 
TDAs. We refer to the different references for details. For illustration, we show here on Fig. [2] 
the TDA A(x,£,t) obtained in Ref. [4] in arbitrary unit; its normalisation would be eventually 
fixed by the experimental data. 

For the purpose of this note, we only need a rough evaluation of the order of magnitude of the 
cross section and will only use the Model 1 of Ref. |4 . When a dedicated experimental analysis 
is being carried out, a careful survey of the cross sections obtained from the different theoretical 
models will be in order. Hence, based on a first experimental study of the £ dependence and 
after having checked the scaling in Q 2 , we shall be in position to see which model describes 
best the physics involved. For this best model, we could then obtain by sum rules relations a 
first measurement of the axial ir° form factor. 

For the following, we shall show results for (ir°\dOAd\ , y) = — 1/2(tt°\uOau\^} expected from 
the different charges of the u and d quarks and using (from isospin arguments) 

(ir+ldOAU^) = (7r°|JC A d| 7 > - (ir°\uO A u\ 7 ). (9) 

This would give Al = l/3A 7r and A£ = -2/3 A" . Note that more realistic models may 
give significantly larger rates. 
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A Q (X,^,t) (in arbitrary unit) 




Figure 2: The 7 — > 7r° axial transition distribution amplitude in Model 1 of Ref . |4 

(for t = —0.5 GeV 2 ) in arbitrary unit. 



3.3 Cross section 

Taking into account the contribution from the fermionic line for the emission by the electron 
of a longitudinal photon, averaging over the real photon polarisation and integrating over ip 
thanks to the 92-independence of the TDA process, we eventually obtain the differential cross 
sectioiE] 



da 



TDA 



dQ 2 dtd£ 9(£+l) 4 <3 8 



W 2 

(-2et)(l-e-(l + 0-^)(Re 2 (^)+Im 2 (X^)). (10) 



For the hypothesis discussed above, the resulting cross section is roughly one sixth of the one 
obtained in [4] for the charged case. The evolution as function of £ is displayed on Fig. [3] Note 
that for small £ (particularly W 2 ^ — * Q 2 ), the cross section shows a peak. 

The Q 2 -behaviour is model independent and thus constitutes a crucial test of the validity 
of our approach. 



A factor 1/4 is missing in Eq.(23) of 0]. 



t=-0.5 GeV 2 s er =20->40 GeV 2 Q 2 =2 GeV 2 

d °V>e,r° n° , 

fb/GeV 




Figure 3: Differential cross sections 



do 



dQ 2 dtd£ 



for the TDA subprocess as a function of £ for 



Q 2 = 2 GeV 2 , t = -0.5 GeV 2 and 3 values of s e7 : 20, 30 and 40 GeV 2 (from bottom to top) 



Conclusion 



We believe that our models for the photon to meson transition distribution amplitudes are 
sufficiently constrained to give reasonable orders of magnitude for the estimated cross sections. 
Cross sections are large enough for quantitative studies to be performed at high luminosity 
e + e~ colliders. After verifying the scaling and the ip independence of the cross section, one 
should be able to measure these new hadronic matrix elements, and thus open a new gate to 
the understanding of the hadronic structure. In particular, we argued here that the study of 
7*7 — > 7r°7r° in the TDA regime could provide with a unique experimental measurement of the 
7T° axial form factor. 
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